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I. INTRODUCTION 



There are a few different starting points for setting of the statistical mechanics problems. 
These starting points are canonical, microcanonical, grand canonical and other ensembles. 
In general, these ensembles are not equivalent. Moreover, nonequivalence of the ensembles 
occurs for most interesting problems, such as new phases nucleations, statistics of small 
systems, metastable states etc. Therefore, statistical mechanics in its present state is not a 
well-posed theory. 

From this standpoint, it is interesting to develop the general methods that can be used 
to different ensembles so as to different models of statistical mechanics. The problem of the 
canonical configuration integral calculation exists in every ensemble of statistical mechanics 
with some variations: for canonical and grand canonical ensembles the configuration integral 
depends on the real temperature, whereas for microcanonical ensemble this integral depends 
on the imaginary temperature The main difficulty in the partition function calculation 
caused the atomic coordinates entanglements in the integrand. Factorization of the integrand 
on the atomic coordinates by means of the functional integral representation solves the 
entanglements problem. But instead of that, the new problem of the non-Gaussian functional 
integral calculation occurs. 

Probably, for the first time the notion of the functional integration was introduced by 
Wiener (2}. In connection with statistical mechanics the functional integration method was 
initialized by Zubarev 0] and Edwards Some mathematical aspects of the functional 
integration method, including approximate calculation of functional integrals, are system- 
atized by Egorov, Zhidkov and Lobanov 

Functional integrals include integration over some set of the auxiliary functions. In ac- 
cordance with domain of definition of the auxiliary functions, there are two basic variants of 
the functional integral representations for classical partition functions in classical statistical 
mechanics. 

1. Functional integration over the functions defined in the real space 0, @, 0]; 

2. Functional integration over the functions defined in the wave- vectors space 0, H, 0, [Io| . 

In both of these variants the partition functions or the configuration integrals calculations 
reduced to the non-Gaussian functional integral evaluations. 

Functional integration approach has essential advantages over the direct analysis of the 
configuration integrals. In particular, this method makes possible 

1. avoidance of the combinatorial problems related to the Mayer expansion 0, @] and 
some other similar problems; 

2. analysis of some simple models statistical mechanics such as the systems with non- 



negative Fourier-transforms of the interatomic potentials [111 . [12|, [13|, [14 



However, the auxiliary field in the functional integrals is an interesting object by itself. It 
gives a new way of the statistical mechanical problems formulation. 

The object of this paper is the formulation of the new variant of the statistical mechanical 
problems setting. This variant based on the functional integration method over the auxiliary 
functions defined in the wave-vectors space. 
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II. FUNCTIONAL INTEGRAL REPRESENTATION IN WAVE- VECTOR SPACE 



The generating functional of a system of iV particles, interacting via two-body poten- 
tial v (r) in presence of an external field <p (r), after integration with respect to momentum 
has the known form 



N , n ^ \ / N 



x exp 



/ * » X 
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, s,s'=l, ■ 
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where A = {2-Kh/mkBT) 1 ^ 2 is de Broglie thermal wavelength, (3 = 1/kgT is reciprocal 
temperature, V is the system volume, D is the space dimensionality, % and k B are the 
Planck and Boltzmann constants, respectively. 

The first exponent in the integrand separates into one-atom multipliers in contrast of the 
second exponent. Suppose that the central interatomic potential v (r) assumes expansion 
into the Fourier series. Then the pairwise interactions energy of the particles in the system 
can be written as 
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where C (k) and S (k) are the collective variables 

N N 



C (k) = cos ( kR *) ; s ( k ) = Yl sin ( kR *) > ( 3 ) 



8=1 



summation over k and k' in (j2j) does not contain the point k = 0, the term — (0) 
compensates the summands with s = s' in the right hand side of (j2j), denote the subsets 
of the wave vector space in which the Fourier transform v (k) of the interatomic potential is 
positive and negative, respectively, 

v + (k) = +v, ke Q + ; ( 
v - (k) = -v, k G IT, 1 J 

w±(k) > 0, n = N/V. 

The second exponent in the integrand of ([!]) as function of the collective variables has 
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the following form: 



/ 
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Applying the Stratonovich-Hubbard transformation 
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for each of the exponents in (jSJ) 
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we obtain the following functional integral representation for Z{<^(r)}: 
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where 



^(^(k), ^(k), {<p(r)}) = Jye 



(V) 



xexp lip ^ J mv+ J^> u (k) [ x + (k) cos (kr) + y + (k) sin (kr)] \ (10) 
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uo (k) is an arbitrary positive function of the wave-vector with the dimensionality of the 
circular frequency, m is an arbitrary positive parameter with the dimensionality of mass, and 
x^ (k) , y (k) are auxiliary variables as a result of Stratonovich- Hubbard transformation. 

Functional T\ (x ± (k) , y ± (k) , {ip (r)}) has a very simple physical interpretation: it is 
the configuration integral for one particle in complex-valued external field ^ (r) : 



* (r) = ip (r) - z 2^ v/ ^7^^ ( k ) [ x+ ( k ) cos ( kr ) + V + ( k ) sin ( kr )] 

" \J mV ~y^ u (k') [oT (k') cos (k'r) + y~ (k') sin (k'r 



(11) 



This complex-valued external field consists of two parts: the "true" external field ip (r) 
and some artificial complex field caused by exclusion of the interatomic interactions and 

represented in the form of the Fourier series with its coefficients \ V^LJ^l. ^ (k) x ± (k), 



mv v ^ uj (k) y ± (k). The expression T\ (x ± (k) , y ± (k) , (r)})^ is a configuration in- 
tegral of a classical ideal gas in the external field $ (r). 

In accordance with (Q, any classical system of interacting particles is statistically 
equivalent to an ideal gas in the complex-valued external field with Gaussian 
distribution of its the Fourier coefficients. 

There are a few variants for functional integral representations like ([9]) of the canonical 
and grand canonical partition functions in statistical mechanics. These representations can 



be used for analysis of some concrete models [ll|, [12|, [13|, [14|]. But, there are not any general 



methods of the functional integrals calculation excepting Gaussian integrals. Therefore, most 
of the results, related to the functional integral representations of the partitions functions 
in statistical mechanics, are obtained by the following ways: 

1. The saddle point method. This method is consistent if Q~ = and leads to Gaussian 
functional integral. The saddle point method relates to the mean field approximation. 
Otherwise (i.e. if f2~ 7^ 0), the location of the saddle point in space of the auxiliary 
field components depends on temperature and Gaussian approximation is inadequate. 
This problem has not solution yet. 

2. The expansions of the functional integral in powers of the auxiliary fields x ± (k), y ± (k) 
with summation of some finite number of the terms in this series. In this connection, 
there is no results on this series convergences. 
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III. HAMILTON! AN OF AUXILIARY FIELD 



The expression (Q can be interpreted as a statistical integral of the system of interacting 
oscillators. Substitute the identity 

, , r /_2 .2 



into (El) and obtain 
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This integral represents the classical partition function of a system described by the Hamil- 



tonian 7i: 
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The multipliers i n the integrand of (fTBl) have the dimensionality of the Planck constant 
and play the role of an elementary cell volume in related phase space. 

Hamiltonian (114p depends on the generalized coordinates x^ (k) , y ± (k) and the general- 
ized momenta p^ (k) , p^ (k) . In addition, Ti. contains the temperature as parameter in the 
term T\ (x* (k) , y ± (k) , Mr)}). 

Thus, the partition functions for a classical many-body system and for the 
auxiliary field with Hamiltonian (1141) are identical. 
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IV. CONCLUSION 



In accordance with relation ffl"3]) . the problem of obtaining the thermodynamic properties 
of many-body system is equivalent to the analogical problem of a system described by 
Hamiltonian (|14|) . This Hamiltonian contains the summands corresponding to a system 
non- interacting classical harmonic oscillators (the first and the second lines in formula ([HP ); 
the third lines of (fT^|) contains interactions between these fictitious oscillators including their 
anharmonicit ies . 

This Hamiltonian leads to some non-linear equation of the auxiliary field motion. The 
analysis of the auxiliary field evolution on base of these equations can be used for the 
statistical properties of the auxiliary fields by analogy with molecular dynamics method. 
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